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ABSTRACT 


/97/i 

It has long been knot’ll that an in^rartant mode of energy loss 
for cosmic ray electrons is i:,v»’'se Compton scatcering with photons 
of starlight. Previous calculations of (dC/dtN due to this process 
have involved non-systematic approximations involving the form of the 
Klein-Nishind formula and the angular distribution of the radiation 
as seen in the electron's rest frame. The present paper considers 
an electron of arbitrary energy in an isotropic therr-rJ. radiation 
field of temperature T. A formally correct expression for <dP./dt)^^ 
is obtained as an asymptotic exi>ansion in the quantity pfeT/CtOgC^)^ 
considered as a small parameter The oft quoted resiilt 
(d^ydt) oc is seen to be the zeroth order term in this expansion. 

/IV 

It is also seen that the energy loss rate c.-anges sign at ; energy 
p, ~ ^ kT as would be expected from thermodynamics. A derivat\'r 
of the zeroth order term is given from classical radiation thecry 
and from this it is seen t ;t thi- term also describes the energy 
loss rate due to synchro tro:\ radiation as well as from inverse 


Compton scattering. 




INTRODUCrriON 


The scattering of energetic electrons hy energy photons, 
called "inverse" tompton scattering, has "been oi as trophy sical 
interest for many years. It was first investigate by Feenberg 
and Primakoff^ as a process by which cosmic ray electrons (an I otons) 
would lose energy during their passage through the galaxy. Later 
Donahue^ applied the general method of Feenberg and Primakoff to the 
case of electrons trapped in orbits about the sun. 

The result of these two papers that the mean energy loss of an 
energetic electron of energy £" is proportional to both the photon 

-A o 

energy density and to c was applied by Hayakawa and Kobayashi 
and by Hayakawa and O '.uda'^ to the problem of the equilibrium of cos- 
mic ray electrons in the galaxy. More recently, Felton and Morrison^ 
have considered this process as a possible soiirce of galactic x-rays®~® 
and gamma -rays®" and Shklovsky^^ has proposed it as a source of x-rays 

in solar flares. 

In the calcuD.ations of Feenberg and Primakoff and Donahue 
the relevant cross section formula is the Klein-Nishina formula 
o(e', X*) for the scattering of a photon of energy ^ ' by a 
stationary electron through an angle X’, In essence the scattering 
probability is expressed in the electron’s rest frame and then trans- 
formed to the laboratory frame to determine the mean energy transferred 
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from the electron to the photon. In the pz'evlous calciilations 

the full Klein-Nishina formula was nv: used hut rather the asyr..ptotic 

forms for e' 4 (Thompson scattering) and for e' -j,. lo. 

Feenherg and Primakoff the two forms are used in the two regions g ' < 
and e* > m^c^ respectively as an approximation to the correct formula. 
Donahue, on the other hand, uses the two formula in the regions e ’ < m^c^/I|. 
and e ' > ^ respectively and connects the two regions of his results 

with an "eyehaJJ ' curve. Both authors assume that the electron is 
energetic enough so that in its rest frame all of the incident radia- 
tion has 6' p=! 0 where rr - is the angle between the photon momentum 
in the electron reso frame and the origineil direction of the electron 
momentum. This assumption obviously limits the validity of the refjults 
to high energy electrons. 

The primary problem with these approximations is that they are 
non-systematic; they do not suggest how to apply a higher order 
correction. In the present calculation the mean rate of energy loss 
- <d£/dt> of an electron in an isotropic radiation field in thermal 
equilibrium is determined. The only approximation used for a wide 
range of is a systematic one in that the resuD.t is obtained as a true 
asymptotic expansion in powers of a small parameter v = C kT/(m^c^)^ 
where £ is the electron energy. and T is the temperature characteristic 
of the radiation field. Wb^ri £ gets so large that z, ~ 1 the expansion 
is no longer useful and tho resul'. must be obtained by methods which 
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are less systematic out which are, nevertheless, quite accurate. We 
will also discover an interesting relationship "between inverse 
Compton scattering, as described by the zeroth order term of the ex- 
pansion, and synchrotron radiation. 
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II FORMULATION OF THE PROBLEM 

Consider an ; 'ectron of energy \ (ve shall express the electron 
energy = y me"', the photon energy s = a mc^, and kT = @ mc^ in 

tenns of the dimensionless parameters ' , o-, and ^ i.ng in a region 

of space in which che photon density, a 0), is gi/cn as a ftuiclion 

of energy a and angle 0 where tf - 0 i.- the pjagle between the photon ar d 
electron velocity •'•ectors. Letting a prime {.') indicate quajitities 
expressed in the electrons rest frame we may write for tha number of 
Compton collisions per unit time 


<dN/dt>^^ = <dN/Ydtn^s. J dn'(G') j aa» (c/y) n* (a», 0‘) at (a*) (l) 

■^o 


where cr‘(a’) is the Klein-Nishina total cross section. Since 
n' (a- ',6) dfi* (S') d<x' is a number density it transforms under the 
Lorentz transformation like an energy ct' so that n* (a', 9) dfi (0')da’/O'* 
is an invariant. 


We may then set 
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n* (a', 0) dfl* (O') dUi* = (a*/ a) n (a,0) dfj (O) da 


and obtain 


J df^(O) J da (a'/Y°-) n (a,0) a' (a*) 


If we denote the lab frsune energy of the scattered photon by a^^ 
■transfer 

the eaergy/ in the scattering process is (aj^ - a) and the mean energy 
loss of the electron is given by 


<“ dy/dt)^^ * ° f (o-'/Y*^) h (o-, 0 ) JdO'(x') o' (a',x' )(ai-cf) (!<.) 
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where c (a',x*) is the differential Klein -Nishina formula for 
scattering a photon of energy oL* through an angle x'- We have the 
usual angle -energy relationship for Compton scattering 


tti* = a’/ [ 1 + a' (l - cos x')l • 


( 5 ) 


Eii 5 )loying the well known Doppler shift formula. 


a’yCl - 0 cos 9x') 

= ai'v (l - P oos 0i') = (6) 

1 + a’ (l - cos x‘) 


(p = - 1 /y = v/c) 


Y (1 - P cos ©i') 

tti - a = a' i 

'' 1 + a'(i 


cos x') 



( 7 ) 
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Consxjlting figure 1 ve obtain the following formula from spherical 
trigonometry 


cos Ox* = cos 0* cos X* + sin 0' sin x* cos 0 . (8) f 

i 

Since the cross section formula cannot depend on 0 we may chose x* \ 

and 0 as our coordinate angles for dfl* (x*, 0) and immediately integrate 
over 0. This has the effect of multiplying equation 4 by 2 tt and re- s 

placing cos 0i* with — q— T where ; 


cos 0x* = cos 0* cos y* . ( 9 ) I 

The Klein-Nishina cross section formula is 



- a - 


ro^ (l + cos^ X*) f 

•(a',x') = — n + a'(i-cos x')l^ ^ 


a'^(l - cos x')^ 

1 + } 

(l+ cos^x’ )[l+c ' (l-cosx' )] 



It shoxild te noted that the variables in this expression are not 
independent since a, a', y, and 0 are related by the Doppler 
shift formula 

a' = Y a (l + 0 oos 0) 
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We now specialize "to the case where n(a^ 6) represents an 
isotropic radiation field i.e., n(a>0) = n(a)ATr 

dn(e) = 2Tid(cos 0) = da'/(0Yo) 

and obtain 


(-(iy/dt^^^ = ^ itro^c j da n (a)/ 


„2 YO-d+P^ p (f^-2f^+2f) ^ 

^ io ^ J (1 ^j'^Tc'^W-ada'a] ;_(-^+2,«f)3 J 


(12) ) 

i 


ipjje integration over f and a* may be done by a strai^tforward and 
repeated application of a good set of integral tables, .\fter some 
time we arrive at the result 



CO N(a) 

<-dY/dt>. =~ Ttro^c r ^ Ha,s) 

•’ o PY^a^ 

F(a,Y) = yC^i^y ) - fi(a/Y)] 

■ a[f2(aY^ ■ ^z(oi/y)’^ 
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where y = + 6^ = y + 

and fi(z) ss (z 4 6 + 5/z) log (l + 2z) 

-(22z^/3 + 2^z^ + l8z + k) (l + 2z) ^ 

-2 + 2 Lia (-2z) ilk) 

fa(z) = (z + 51/6 + 5/z + 3/2z2) ^og (i + 2z) 

- (22z^/5 + 28 z^ + 103z/3 + 1? + 3/z) (l + 2z)"^ 

- 2 + Lia (-2z) . (15) 

!Ihe function Iii 2 (z) is the Eulerian Dilogarithm^^ defined hy 
z log(l - z’) 

LigCz) = - dz* for coniplex z 

o 

00 

“I 

nsl 


z] < 1. 
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III EVALUATION OF THE INTEGRAL 


We vill now conjider that the photon deuoity is described 


reasonably well by the RLaock radiation fonnula 






where (a) * J ® ^ “ "energy” density and ^(p) is the Riemann 


Zeta function defined by 


= 1 + 1 / 2 ^ + 1 / 3 ^ + . . 


and c (^) = ttV 90 


Expression 15 may now be written 


<dy/dt> = nro^c <a> N (y) 


(15/rf*) ri » fi(z) dz 

N(y) = r z I y f 

0‘*2(v^-l)'sU **0 exp (z/^0) - 1 **, 


fx(z) d2 

o exp (zy/©) - 1 
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(15/tt^) , 

rL 1 

, zfa(z) dz 

- - - — ^ 

" zf 2 (z) dz 


J 

• Y 

0 exp (z/^©) -1 

0 exp (zy/®) - 1-* 


(IT) 


where we have transformed varjahles in such a way as to make the 
arguments of fi and fa the varr.ahle of integration. 

Up to this point our calctilations have heen exact within the 
firework of the physical situation that we have considered. We 
now ask whether there is some systematic approximation scheme 
that will allow us to evaluate the integrals in equation I7 in some 
siniple manner. 

To this end we notice that f^ and fa have a second order pole 
and a branch point at z = - I/2 (a remnant of the pole in 
a(a',x') at cos x' = 1 + l/2')* This '»ieans that the power expan- 

00 CO 

sions of fi and ^1 = are convergent 


n=d. 


n=a 


only for |z| < I/2 and our integrals aie over the range 0 a z < » . 
However, we notice that the term [exp (kz) - 1] ^ is a 

function that peaks at l/k and drops off as exp (-kz) for values 
of z significantly greater than l/k. Therefore if ©7- < I/2 
only the portion of fj_ and fg for z < I/2 will contribute signifi- 
cantly to the integral. Since © = kT/m^c^ « 10 ® for T = 6,000°k 
we may consider y & as a small parameter of order t (since y s 1, 
©/y ^ y 0) . The expansion coefficients and for f x 
respectively are obtained, in a straLL^tforward manner from the 
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knovm expansiors of log (l + 2z), (l + 2z) and Ll 2 (- 2 z). These 
coefficients for n up to ten are given in table I. 

Zf we now insert the series form of f^ and fa into the inte- 
grals in equation 17 and ignore the fact that they are not con- 
vergent for jz] > 1/2 we will obtain a formal series for NC-y) 
which we hope will not be in error by very much. Making use of 
the formula 


r 


0® Uj 

X dx 

exp U/c)-l 


= n' C (n+1) 


n+1 

e 


( 18 ) 


we have 


“ fi(z) dz 

f ■ , , ' 

j exp (z/e) -1 


I 


A^n! C (n+1) 


n+1 

e 


00 , V 

z fa(z; dz 
o exp(z/e) -1 


= Vb (n+l) ; Q (n+2) 
L. n 
n 


(19) 


We may see at once that these* series are not convergent since from 
the known circle of convergence of the series for f ^ and fg 



lim A 


we have 


n 


n-'>® A. 


h-1 


lim B 

_r 

n-<B B 


= 2 . 


n-1 


lim 

therefore since co C(z) = 1 we have 
z-» ^ 


n+j. 


lin ^(n+1)^ 

n-’" (a Jn-1): C(n)eV 


= 2n« • 


So for any finite g there is a value of n -- 1/2 e beyond which 
the terms of the series grow without limit. It is fairly easy to 
see that this is the same series that would be generated by 
repeated partial integrations of the j.ntegrals in 19 since 
A^n.' = d^fi/dz^. Furthermore, it is demonstrated in the ap- 

pendix that this is in fact a correct asymptotic expansion of the 
integrals in the sense that 

N 

F(e) = + 

n=l 

andRjj(e) = 0 . 

Noting from bable I that 


Ao s» Aj_ SB Ag » Bo *= Bj_ SB 0 
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we may, after some rearrangement of terms, write for W(y) 


K(v) = (15V) [v"a^ 5 C_^ 5 (y) 

m=o 


■®m+2W C (dh-^) (y®)“ 

where 

“ 1-1 

®n<^> 'I G ) (1 - y'^) = 2y“-^(y^-D* 

issodd 


¥e see that for «^© ^ 1 (in our case y ^ 10®) this series gives 
an excellent approximation provided you don’t sum heyond n tw 1/2y®. 
The zeroth order term in this series is 


(y^ A® C® (y) - B 2 Ca (y)) 5.’ 


C (4) (25/tt^) = (A (V^-1) ‘2Ba) 

3 


= H-Aq (y^-l) 


since 2Bg = 3As. 




k 
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Inserting this approximation to N(y) in expression 16 we have 
to zeroth orier 

= 5.555 TTTo^c p(p/m^c)^ (21) 

where p is the photon energy density in conventional units. We 
see that for electrons of sufficiently low energy such that pc <f 
the energy loss rate is proportional, to p^ rather than£*^. 

We note further that expression 20 is negative for y sufficiently 
close to 1. This means that a very low energy electron car, on the 
average, gain energy from the radiation field. Setting expression 20 
to zero we have to first order 

A34.' C W = - [^A4 - 3B3 + (y^-1) (8A4 - 

X kl C (5) (y©) (22) 

Since equation 22 states that y^-1 = 0(e) we may neglect the term in 
y^-1 on the right hand side and obtain 

(3B3-4A4) C(5) 

cW 

« 5y® (C(5)/C(^)) 


0.958 (5y®) 


(25) 
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We see that this is a few percent below the equipartition energy 
for a relati'.'istic gas^* Y^'l = 3 yS- A similsu: calculation for the 
energy loss of a test particle in a Maxwellian, hard-sphere gas 
gives {d(f/dt\ = o for = O.98I (5kT/2 ) which is in qualitative 
agi’eecent with equation 25- 

For ^ 'tile eisynrptoti^ series 20 is no longer useful. In 
this region we mast resort to less systematic hut nevertheless 
quite aceorate methods of approximation. First of all note that if 
y3 ^ 1 then y/S ^ 1. This means that the second and fourth integrals 
in expression I 7 may he completely ignored. We may, in fact, rewrite 
expression I 7 , noting that y 2y, as 


. . fi(z) dz . - t» z fo (s) dz -1 

»(v) = (15/n ) [(H ) exp(z/2Y®)-l ■ > f S5 (z/2v0) - 1] 


( 21 .) 


Prom expressions 14 and I 5 we also see that hoth f ^ and tz tend 
to z log (az) for z large where a = 2e in fact for 

z s z' = 7*75 X 10* this approximation is good to within one part in 
10*. We may, therefore, perform the integrations indicated in 
expression 2k in two parts; from zero to z* we evaluate the inte- 
grals numerically and from z' to * we use z log (eiz) for fi and fg 
and obtain analytic expressions. 

Carrying out this procedure we obtain the following, rather xm- 
transparent expression where T - 2€>y 


N(y) = (6o/tt*) 


z' 



fi (z) dz 

exp(z/r)- X/ numerical 
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- r ^ z' log (az') log (l-e 


+ r’^ (log (az') + 1) Lis (e“^/^) 


+ r~“ Y -El(-nz'/r) 1 


nsl 


.z' 


1/2 [ r'^ ( r z fa (z) dz ^ 

' 0 exp (z/r) -1 aumerical 


- r® z'^ log (ii^') log (1 - exp (-z'/r)) 


+ r z' (2 log (az') + 1) Ida (exp (-z*/r )) 


+ r ^ (2 log (ax') + 5) Ida (ej^ (-z'/r)) 


4- \ n't 

L n3 J J 


n=d 


(25) 


Id 3 (z) is the trilogarithm where in general 
z Li„ ,(z') 

Li (z) 5= r dz* 

n J » • 


« ) z /m 

msd 


Z < 1 


I 
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and Ei(z) is the exponential integral defined by 


- Ei(z) 



dt . 


The values of this expression mqy be calculated quite easily on 
a coniputer; however, it is instructive to consider the situation 
vhen Y (hence r) becotsss very 2arge. A brief examination shows 
that every term in the expression tends towards zero with the 


exception of the first series ) — — — 

properties of Ei(-x) we have 


From the known 


Ei(-x) = C + log (x) + f(x) 


where C = 0.577215669 is Euler's constant and f(o) = o and 
|f(x)| jC + log(x) ! as X 


V f ^ nx) 

Therefore ■ converges uniformly for all finite x and 

converges to zero for x equal to aero. We are then left with 


C + log ^ 

-I — -y. 

n n2 


n 


= C ( 2 ) log 




+ const. 
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We then have 


N(y) - (l5/'n‘^0^) (C(2) log V + const.) 


(26) 


as y becomes very large. 
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IV RESULTS AND CONCLUSION 

In figure 2 we present curves of N(y) = - (d£/dt^(TTro^o)”^ 
as a function of y for a range of T from 5>000°K to 

10,000°K. The necessary computation was done on an IBM 709^* 

For values of 2y0 ^ 10 ^ the asymptotic series 20 was used in- 
cluding the seventh order term. For values of 2y0 > 10 ^ ex- 
pression 25 was evaluated. 

It is immediately seen that the approximation - (d£!/dt) 
is very good for 4 < y < 4,000. Below 4 the dependence on p^ 
rather than^^ is manifest. Bet\.een about 10® and IQ® the curves 
achieve the form of expression 26 namely proportional to log^* 
and T The values of T'Cy) for y -• 1 can not be shown in figure 
2 due to the logarithi-ic scale, however, the zeroth and first order 
terms of expression 20 give q\iite accurate values in the region 

y ~ 1. 

It is Of interest to conipetre inverse Compton scattering 
and synchrotron radiation as an energy loss mechanicsm for ccomic 
ray electrons. To this end we first note that for small values 
of y 0 the inverse Compton scattering process is a classical radia- 
tion process; in this limit the Klein -Nishina cross section is just 
the Thompson scattering cross section. 

For an electron in arbitrary electromagnetic field the inst€ui- 

15 

taneous radiated power is given by 



I t 
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2 e" 4+ 8 lX H)2 - 

p = J 

3 c® ^ 1 - 

4 

2 ® ^ (e X E)2 + (o X H)2 - 2« . (E X ^ 

= (e 2 , ± — r_ri 

3 ^ 1 - (27) 


If we now assume that any energy flow is isotropic ((p.(ExH))=0) ? 

A.' I 

and that the fields are unpolarized (<(B x E)^) = £ E^; ' 

^ >»v 3 ? 

((b ^ H)^) = £ b^ H^). ¥e now have 


<P>^^ = I Ttro^ (cE 74 tt) + 3.555... nro7 i-} (28) 

8tt me' 


If we consider the situation that the only ^ fields present are 
radiation fields the radiation energy flux incident on the electron 
is ^ixst ® TtTo^ (cE^/4tt) so that the loss of mechanical energy is 

3 

Just radiation out minus radiation in or 

- <d£/dt>^^ = 3.555 ... T, ro^ C ( -g;; ) ( ^ ) (29) 

idiere the energy density includes reidiation and static magnetic fields. 

This is, of course identical to equation 21 and we see that the 
relative importance of inverse Compton scattering compared to synchro- 
tron radiation depends only on the energy density of the radiation 
field versus the energy density of the magnetic fields^. We see, 
therefore, that in the galaxy where the energy density of "both starDJLght 
and magnetic fields is of the order of 1 ev/cc the two processes will he 
oh a rou^ly equal footing. 
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APPENDIX 

GENERATION OF AN ASYMPTOTIC EXPANSION BY REPEATED PARTIAL 

INTBG . iTIONS 


Consider the integral 

^ " j exp(z/e)-l 


vhere f(z) is smalytic on the positive real line including zero 

and it and all of its derivatives increase no faster than a poly- 

“ / 

nomial as z We may expand [exp (z/e;) -l^ as > e idiich 

n^ 

series converges uniformly in z for z >■ o. Since the series does 
not converge for z = o we must also demand that f(z) go to zero 
at least as fast as z for z -• o. Since the integral now has con- 
tributions only in the region of uniform convergence we may inte- 
grate terra by term. Integrating a particular term by parts N 
times we obtain 



-nz/fi 


N 


dz 




m+l 


m=o 


where R^^ = 


r“ d:! 


Due to the analyticity of f(z) on the positive real axis and the 
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Htni t.ftt.Ton on its growth as z -• co f nay be bounded by 


|f^^"’^^z)l ^ Ajj t so that 


. . N+2 J 


Nt2i-q. 


If we now stun over n remenibering that f(o) = o we have 


N 




Ift=l 


where = C f^™^(o) 

and iR’jjj s Ajj C(Nt2) ^ (U+2tq) e 

IT 

This expansion, therefore, satisfies the definition of an 
asymptotic expansion of I as e -* o* 


I I 



n. 
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TABLE I 

Expansion Coefficients for f^ and f^ 


n 

n 


0 

0.0 

0.0 

1 

0,0 

0.0 

2 

0.0 

0. 13333320 x 10 

3 

0. 88888827 x 10° 

-0. 28444412 x 10^ 

4 

-0.27999992 xlO^ 

0. 69999923 x 10^ 

5 

0.78400015 xlO^ 

-0. 16784715 X 10^ 

! 6 

-0.20520631 xlO^ 

0. 39161810 X 10^ 

7 

0.51156461 xlO^ 

-<). 89469178 x 10^ 

8 

-0.12304761 xlO^ 

0. 20112988 X 10^ 

9 

0.28815807 xlO^ 

-0. 44645477 x 10^ 

10 

-0.66115232 xlO^ 

0. 98096851 x 10 







Figure Captions 


Figure 1. Ar.gles involved in the scattering process viewed rrom the 
electron center -of -mass system. 


F'gure 2- N(«y) versus y for values of © corresponding to T = 5j^'00^K, 

7,0C0% and 10,000^. 
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